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Introduction 

In the L = 2 to 5 range in the outer radiation belt, Davis, 
Hoffman, and Williamson [ 1964 ] find that the spectra of the 
relatively stable 0.1 to 5 Mev protons show smooth but large variations 
with L and equatorial pitch angle, or , as in Fig. 1. Protons near 
the earth and at a Q near 90 ° are more energetic than those at 
larger L and at smaller ot Q . The spectra are well represented by 
e‘ E/E °. 

Inspection of the experimental data showed that E 0 « L “ 3 « B. 

This was very suggestive since to have a characteristic energy of 
the system vary with the magnetic field suggested betatron 
acceleration and suggested that the particles were moving radially 
across the field. Such a process must violate at least the third 
adiabatic invariant of motion since the particle moving radially 
does not conserve the magnetic flux inside its orbit. 

Kellogg [1959] first suggested that the radiation belt might be 
formed through magnetic disturbances in which the third adiabatic 
invariant of trapped particles is violated without violating the 
first and second invariants. Violation of the third invariant 

* The first part of this paper follows closely the material in 
Nakada, Dungey and Hess, J.G.R., JO. 
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allows motion in L-space. As particles move closer to the earth 
obey tend to gain energy with the maintenance of the first 
invariant since, for example, E/B is a constant for 90° pitch 
angles. So, this process can introduce acceleration of protons. 
Kellogg 1 s suggestion has teen adopted for this study although the 
mechanism for motion in L space is unspecified. It has further 
been assumed that motion in L-space is rapid compared with atmospheric 
loss and scattering processes except very near the earth and that the 
geomagnetic field is sufficiently well represented by a dipole. 

Energy and Angle Variations 

If the first and second adiabatic invariants of trapped particles 
are maintained during motion in L space, changes in both the energy 
and equatorial pitch angle can be calculated. The first invariant 

is 


P 


E sin 2 q n _ E L 3 sin 2 gp 
B 0 “ .312 


( 1 ) 


where Bq is the equatorial magnetic field. The second invariant is 


- m 


<f> v cos a d£ a 


j 


mv cos a 1 - JH L F(or c ) 


( 2 ) 
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where m is the mass, v the velocity, o' the local pitch angle and 
l is along the guiding center. The integration is over a complete 
north-south oscillation. 

Since p and J are constants, equation (l) may be divided by 
me square of equation (3) to give 

L t " fXg'q? 1 = constant - (3) 

From this, the changes in ot Q with L have been found and are shown 
in Figure 2. As Davis and Chang [1962] have indicated, particles 
diffusing inwards assume flatter helices and move closer to the 
equator. 

These changes in ar 0 with L and equation (l) may be used to 
find the variation in energy with L and ot 0 . Results are shown in 
Figure 3 for protons having ao values at L = 7 as indicated on the 
curves. Energies are relative to energies at L " 7. 

The energy transformation is given by 

E L 3 sin 2 a = E L 3 sin 2 or ( 4 ) 

oL s s os ' ' 

where the subscript s refers to the source. If the source spectrum 
has an exponential form, the transformation is given by 
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e -E s /E 0S 


- E L 3 sin 2 gn L 
e E os L s 3 sm 2 a os = e “ 



(5) 


which shows 


E 0 L 3 sin 2 Qf r = E L 3 s i n 2 a 

u OL OS S OS 


( 6 ) 


Fran this it can be seen that an exponential source remains 
exponential after L space motion and that Eo changes in the same 
way with L and a 0 as has been calculated for a single particle in 
the previous section. 

These two predictions of the model may be compared with 
experiment. The first prediction, that the spectra retains its 
exponential form, is in agreement with experiment. To test the 
second prediction, measured E 0 [Davis et al., 1964] have been 
plotted in Figure 4 as a function of L with appropriate changes in 
a 0 with L. The labels on the curves refer to <*o values at L = 7 • 

The dashed curves in Figure 4 are taken from Figure 2 for corresponding 
changes in E with L and ar 0 . The changes in Eo with L show good 
agreement between the model and experimental results. The experimental 
results also give the same trend as the model in the change in the 
slopes of the curves with u 0 . 
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If the dashed curves in Figure 4 are extended, they intersect 
near L = 10. This intersection is where the spectrum is independent 
of a 0 and thus gives a source location with the simplest assumptions 
about the source. 

Flux Variation with L 

Now let us consider the way the flux of particles will vary with 
motion in L. The density of particles in phase space, f, is given by 


f(v,u) 


dN 

dA dt p 2 dp dQ 


( 7 ) 


But, the pitch angle distribution in protons/cm 2 -sec-ster-Mev measured 
by Davis, Hoffman and Williamson is 


j(E,of) 


dN 

dA dt dE dO 


( 8 ) 


This gives 


f( v a ) = j- (E ^ cy) ^ 

uv,ar; p2 y dp 


( 9 ) 



This shows the simple relationship between the measured fluxes and the 
velocity distribution function. Our model of the drift process has p 
and J constants of the motion. We can study the way the particle flux 
varies with position by studying 


f (l- 1 !) ) > Oii (p^* 3 


The values of f were computed as functions of L for the data of 
Davis and Williamson for many pairs of the values of p and J, and Figure 
5 shows f plotted against L for fixed p and J. The curves for all p and 
J were normalized in the region L = 4 to compare their shapes and the 
following interesting empirical fact shows up 


f(p,J,L) = g(L) h(p,j) 
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That is to say, the function f is separable and the L dependence is 
essentially the same for all p and I. This being the case we can 
write 


f(p n ,L) 

f( M-i j 



C(L) 


( 12 ) 


where C(L) can be read directly from Fig. 5 and it can be used to transform 
.xes from one L to another. -Fig. 5 shows f varying considerably with 
L. It is seen that Of/QL) j is always positive, suggesting that the 
pur side source is at large L, the particles duffusing inwards and loss 
processes reducing f further in. The small slopes of Figure 5 at the 
larger L values implies that loss processes are probably relatively 
un hap or "cant there, the slope probably being due to diffusion of particles 
away from the source at the outer boundary. The much larger slopes 
at the lower L values imply that loss processes are relatively important 
in this region. 

We now have a scheme for transforming proton fluxes from one position 
in 3,L space to another position. Using equations (l) and (2) we 
know how the energy E 0 of an exponential distribution varies and how 
the equatorial pitch angle changes. Using the data in Fig. 5 we bave 
an empirically determined scheme for transforming particle fluxes as 
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given by equation (12) . Starting with the data of Davis, Hoffman and 
Williamson shown in Fig. 6 for equatorial particles at L = 5> we 
have transformed this to L = 2.4. This equatorial data has been fit 
by an expression 


J ± (E) 


= 2 a. e 

i 
x 


-E/E 


oi 


_ 7 -E/.088 - = 

= 3.8 x 1.0 7 e ' + 1.65 x 10 5 


r E/ * 465 (13) 


Obviously the second term in equation (13) is not very well known since 
from Fig. 6 there are only two 'points that can be used to fit this. 

This transformed data is compared with the Relay 1 data of Fillius and 
Mcllwain in Fig. 7* Over the range where comparison is possible the 
agreement is very good. The transformed flux and the experimental flux 
an L □ 2.4 agree to within a factor of 2 which is as accurately as 
CU) is known from Fig. 5- Beyond E = l8 Mev the comparison is impossible. 
Davis' data runs out here. The experimentally measured fluxes at 
L = 2.4 fall significantly above the extrapolation of Davis' data shown 
dotted in Fig. 7 suggesting that a third term should be added to the 
flux description in equation (13) • For part of the following discussion 
we have added a third exponential 


-E/4.95 


(14) 


J 3 (E) = 1.0 x 10 5 e 
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Off-Equator Protons 

Having found reasonable agreement between measured fluxes and 
transformed fluxes for equatorial particles, let us now examine off- 
equator particles having ao ^ 90°. We start with a fit to the Davis, 
Hofftoan and Williamson data at L = 5 similar to that given in 
equation (13) but with values of the constants a^(a 0 ) and E Q ^(a 0 ) 
varying with the pitch angle. Davis has made such a fit to his data so 
this initial data is well known. Using the, same transformation as before 
involving equations (l), (2) and (12) the flux measured at L = 5 has 
been transformed to other L values and is compared with the Relay 1 
data of Fillius and Mcllwain for off -equator locations in Figures 8 to 
12. In these figures the transformed flux curves are shown labeled by 
which exponential 1, 2 or 3 from equations (13) or (ll) dominates in 
determining the flux above a certain energy at a certain location. 

From inspecting Figures 8 to 12 we can arrive at the following general 
statements : 

(l) Wear the equator the agreement of the measured fluxes 
and transformed fluxes is quite good where direct comparison can be made. 
This agreement remains quite good going off -equator up to a location 
dependent on -energy and location but given very roughly by B/B 0 ~ 3 
for proton energies of a few Mev. For higher energies (Fig. 11 and 
Fig. 12) the region of agreement is quite small. 
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(2) For locations well off-equator there is a decided 
discrepancy between the measured fluxes and transformed fluxes. The 
measured fluxes are consistently higher. 

There are at least two ways in which the off-equator fluxes might 
be made larger than those given by the transformation here. First 
there may be processes which violate either the first or second 
adiabatic invariant that usually will tend to move particles down 
field lines . In this way we could populate off -equator locations by 
moving particles from the equatorial region down field lines. We 
must be careful not to do this so efficiently that the equatorial 
fluxes are changed significantly, otherwise we would lost the agreement 
chat we do have at the equator. But most of the Davis protons in a 
cube of force are confined quite close to the equator so maybe supplying 
die lower altitude population can be accomplished without damaging 
the equatorial agreement. 

A second way of supplying the off -equator proton flux, which seems 
more likely to me, would involve a second source of protons. It is 
well established that neutron decay protons produce most of the 
observed E > 30 Mev protons for L < 1 . 5 . The neutrons decay source 
extends outwards into the outer belt, falling off about a R" 2 . These 
protons will also be acted on by the process that violates the third 
invariant and these protons will also move in L. These will tend to 
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diffuse outwards away from the source. Tverskoy (196k) suggested that 
vhis process might be responsible for the sharp outer edge of the 
inner belt proton distribution. We have the interesting possibility 
that outer zone protons diffuse inwards to fill the region near the 
equator at low L values and simultaneously inner zone protons are diffusing 
outwards to fill the higher B parts of field lines as shown in Fig. 13 
This is what is expected by analogy to gradient diffusion where particles 
tend to move away from regions of high flux. This idea has not been 
tested quantitatively. 

Other Particles 

Recently experiments have indicated that besides protons, other 
particles also undergo L diffusion. Frank (1965) has shown a radially 
inwards moving wave of electrons in the outer belt following a magnetic 
storm. The wave had a radial velocity 


v as kL 8 


( 15 ) 


This wave has the right properties to be due to L diffusion concerning 
and J. The magnetic pumping process described bv Parker (19SO) 

Davis and Chang (1962), and Nakada and Meaa (1965) will operate equally 
well on all particles that have the same drift period around the 
earth since this is the only particle diameter that enters the theory. 
This means that electrons and protons of the same energy will experience 
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roughly the same L diffusion. However the electron problem appears 
i*.cre complicated than the proton one. Other processes such as precipitation 
or pitch angle scattering and in general short lifetime seem characteristic 
of electrons, while for protons L diffusion appears to be dominant. 

Van Allen ( 1965 ) has recently identified a-particles in the 
outer radiation belt. If some of the protons in the outer belt are due 
to L diffusion inwards from the magnetopause, then we would expect 
a-particles there too. It is now quite well established that there are 
several percent a-particles in the solar wind. If these get reasonably 
therraalized in the transition zone, as seems to occur for protons and 
electrons, then we should have for the flux of a-particles in the 
transition zone 


j^(E) & k j (E) where k ~ .05 (16) 


Assuming that these a-particles can get through the magnetopause, 
as seems to be the case for the protons, then by analogy to the known 
proton energy spectrum inside 


j (E) = k e" E p /E °P 


(IT) 
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ve v/oulu expect the a-partiele energy spectrum to be 


Of 


(E) 


k e -20 '^ 00, 
a 


(18) 


Also as a result of the thermalization we should have 


TP - TP 

o a op 


(19) 


for all locations in the magnetosphere where loss processes ar.e not too 
important. Loss processes will be different for protons and a-particles 
and probably will change the two spectra in different ways. Loss 
processes seem important for protons for L < 3. 

This situation where the proton and alpha particle energy spectra 
are the same should be true as long as we have steady state, without 
loss for both particles, that is, where the Fokker Planck equation can 
be written (Davis and Chang, 1962) 


_d 

dr 


[<A r > f] 


+ * 


3r2 


df 

at 


[<( A r ) 2 > f! 


0 


( 20 ) 



‘fhe spectra should he the same for steady state because we would 
expect the diffusion process to operate on the two particles in similar 
ways so that 


<Ar > <Ar > 

a - P 

<( A r ) 2 > " <( Ar ) 2 > 

a p 


which from (20) would give j (E.) = k j (E) . The two sets of particles 

"a p 

don't have to diffuse with the same velocity. They can move through 
each other and probably will,, because protons and u-particles of the 
same energy will have different drift periods 


t a 



r me 
L zeBR 

e 



v/) ! 


-1 


From (22) for protons and a-particles of the same kinetic energy 
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■Inis me an s that a different set of magnetic disturbances will cause 
violation of the § invariant for protons and a-particles of the same 
energy. 

A prediction of this model of L diffusion is that inside the 
magnetosphere the proton and a'-particle energy spectra will be related 


4 (E> 

V® 


k e 
a 


C/Ec 


k e 
P 


"E/Eo 


k ~ .05 


( 24 ) 


in the region of L space where loss is unimportant (roughly L > 3). 

Also if the suggestion put forward here is correct, that outer belt 
protons at large B/B 0 , i.e. well off the equator, are drifting outwards 
from the inner belt, then we would not expect a-particles to exist in 
this region of space. If a-particles are found at large B/B 0 , then 
maybe violation of the p or J invariants is important for outer belt 


protons . 
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FIGURE CAPTIONS 


Lov energy proton spectra at various distances from the earth, 
fro.;. Davis and Williams on, (1963^ 

Variation of equatorial pitch angle with L for motion conserving 


5 . The ratio of proton energies at L to that at L - 7, assuming p and 
J constant, for various equatorial pitch angles at I = 7 • 

k. Variation of E 0 with L for various equatorial pitch angles. 

5- The velocity distribution function f vs. L for various values of 
P and J determined from Davis data. The curves for different p 
and J have been normalized in the Region L ~ 4. 

The exponential fit to the energy spectrum of Davis, Hoffman and 
..■'ill iamson, 1964. 

X- comparison of the transformed Davis, Hoffman and Williamson data 
with the spectrum measured by Fillius and Mcllwain at L = 2.1. 

6 . A comparison of transformed Davis, Hoffman and Williams on ,^1961^ 
data with measurements by Fillius and Mcllwain of protons of 


1.1 no 14 MeV. The curves arc identified by which exponential 
term from Eq. 13 or 14 dominates in the transformed data. 

A comparison of transformed Davis, Hoffman and Williamson, (196!^ 
data with measurements by Fillius and Mcllwain of protons of 1.6 
to 7.1 MeV. The curves are identified by which exponential term 
from Eq. 13 or 14 dominates in the transformed data. 



A comparison of transformed Davis, ID'.ffma n and Williamson, \ v I$io^ 
data with measurements by Tillius and Mcllwain of protons of 
2.25 to 1.7 MeV. The curves are identified by which exponential 
ten;; from Eq. 13 or lb dominates in the transformed data. 

A comparison of transformed Davis, Hoffman and Williamson, Cl9^*-y 
data with measurements by Tillius and me II wain of protons of 
18.2 to 25 MeV. The curves are identified by which exponential 
term from Eq. 13 or l^i dominates in the transformed data. 

A comparison of transformed Davis, Hoffman and Williamson, i 19^' 
data with measurements by Tillius and Mcllwain of protons of 
35 to 63 MeV. The curves are identified by which exponential 
ter;.- from Eq. 13 or lb dominates in the transformed data. 

A possible mcdel to explain off -equatorial low energy protons. 
Protons travel in from the magnetopause near the equator and travel 
outward from the inner radiation zone off -equator. 
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